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Abstract. There are thousands of confirmed detections of star forming 
galaxies at high redshift (z > 4). These observations rely primarily on 
the detection of the spectral Lyman Break and the Lymana emission line. 
Theoretical modelling of these sources helps to interpret the observations 
in the framework of the standard cosmological paradigm. We present 
results from the High-z MareNostrum Project, aimed at constructing a 
panchromatic picture of the high redshift galaxy evolution that will im- 
prove our understanding of young star forming galaxies. Our simulation 
successfully reproduces the observational constraints from Lyman Break 
Galaxies and Lymana emitters at 5 < z < 7 . Based on this model we 
make predictions on the expected Far Infrared (FIR) emission that should 
be observed for LAEs. These predictions will help to settle down the ques- 
tion on the dust content of massive high-z galaxies, an issue that will be 
feasible to probe observationally with the Atacama Large Millimetre Ar- 
ray (ALMA). 



1. Introduction 

Observationally, there are two main ways to detect high redshift star forming 
galaxies. The first technique uses broad band measurements in order to detect the 
drop in flux due to the Lyman Break in the galaxy's Spectral Energy Distribution 
(SED). Galaxies detected through this technique receive the name of Lyman 
Break Galaxies (LBGs). The second technique aims at detecting the Lymana 
emission line either by narrow band filtering techniques or by direct spectral 
measurements. Galaxies detected through this selection technique receive the 
name of Lymana Emitters (LAEs). For some of the high-z LBGs there has 
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also been detection in their rest-frame near-IR which has allowed to provide an 
estimation of their stellar masses. 

In this work we describe our efforts to build a panchromatic picture of 
high-z galaxies at 5 < z < 7. Our work is based on the results from a large 
hydrodynamics simulation that follows the gravitational collapse of dark matter 
and gas and including radiative processes and the ensuing star formation, to- 
gether with associated chemical and mechanical feedback from exploding stars. 
Based on the analysis of the numerical galaxies, we have suggested a model to 
account for the dust extinction and the radiative transfer of the Lymana line. 
We have shown how this model is in fair agreement with the observed properties 
of LBGs and LAEs (Forero-Romero et al. 2010,2011). Special attention was 
made to reproduce the observed fraction of LBGs with a strong Lymana emis- 
sion (Forero-Romero et al. 2012). In what follows, we will describe the main 
features of our model and will conclude by pointing out the expected trends for 
the FIR emission that will be probed by ALMA in the near future. 



2. The Simulation and Galaxy Finding 

The High-Z MareNostrum Simulation follows the nonlinear evolution of a cubic 
region of 50 hT 1 Mpc (comoving) on a side. The dark matter and gas distri- 
butions are sampled with 1024 3 particles for each component. The gas physics 
includes radiative cooling and photoionisation from an homogeneous UV back- 
ground switched on at z = 6. Star formation and feedback (chemical and dy- 
namical) are included following the Springel and Hernquist (2003) model. We 
identified all gravitationally bound objects using the Amiga Halo Finder (Knoll- 
man &; Knebe 2009). All objects with at least 1000 particles are kept for further 
analysis. More details can be found in Forero-Romero et al. (2010). 



3. Modelling the UV Continuum and the Lymana line 

The stellar SED is constructed from the stellar particles contained in the galaxies 
detected as described in Section [2TJ The main assumption under this calculation 
is that each stellar particle in the simulation can be treated as a burst of star 
formation with a given mass, age and metallicity. We assume a constant Salpeter 
IMF. 

We implement a phenomenological model to quantify the extinction pro- 
duced by dust in the stellar SED. The physical model assumes that extinction is 
divided into two contributions that affect different stellar populations of different 
ages. Old stellar populations see an effective extinction from an homogeneous 
Interstellar Medium (ISM), while young stellar populations suffer additional ex- 
tinction due to the molecular clouds where they are embedded at birth. 

The calculation of the escape fraction of Lymana radiation assumes the same 
geometry as the UV continuum extinction. We have developed the Monte-Carlo 
radiative transfer code CLARA to follow the path of Lymana photons through 
the ISM in simplified geometrical configurations. Details of the implementation 
of this model can be found elsewhere (Forero-Romero et al. 2011). 



Panchromatic view of high-z galaxies 
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4. Results 

• LBGs, LAEs and Stellar Masses: Our numerical results for the lumi- 
nosity functions for LBGs are in good agreement with the observational 
results. Only a minor disagreement is noticeable at the faint end of the 
luminosity function, where there seems to be an overabundance of simu- 
lated dwarfs with respect to the observations. UV colours as parametrised 
by the (3 slope as a function of restframe UV magnitudes are also in agree- 
ment with results derived from observations. This shows that the dust 
correction we have applied to the simulated galaxies is reasonable (Forero- 
Romero et al. 2010). The simulated LF of LAEs shows a good agreement 
with observations at the bright end, while showing an overabundance at 
the faint end, which is somewhat more noticeable than in the restframe 
UV (Forero- Romero et al. 2011). Recently, we have shown that the stellar 
masses as a function of intrinsic UV luminosity are also in good agreement 
with the observational estimates, at least for the most massive systems 
(Forero- Romero et al. 2012b). 

• LBGs as LAEs: As part of the international observational effort to probe 
the reionisation epoch using LAEs and LBGs, a new kind of observations 
were conducted. Galaxies detected primarily as LBGs were then followed 
up spectroscopically. From these measurements one can determine whether 
a LBG is also a LAE by requiring the equivalent width to be larger than 
a fixed threshold value. Using these results we can consider now a new 
kind of statistics, the fraction of LBGs showing strong Lymana emission, 
AL yQ . The results from our simulations, without any further change, fully 
agree with the observational measurements of this fraction. We have also 
shown that the Lymana escape fraction is decreasing with galaxy mass 
(Forero- Romero et al. 2012). 

• FIR Emission For each galaxy in the simulation we calculate the ex- 
pected flux measured in the ALMA band at 353 GHz as a function of ob- 
served Lymana luminosity. We find that the brightest LAEs do not have 
a high FIR luminosity, passing completely undetected by ALMA (Fig. 1). 
This prediction of our model can be quantitatively understood as follows. 
Bright galaxies in the FIR have a high star formation rate and dust mass 
contents. These high dust values and neutral hydrogen associated with the 
most massive FIR bright galaxies naturally corresponds to very low escape 
fractions. 



5. Conclusions and Outlook 

We have constructed a model for high redshift galaxies (LBGs and LAEs) based 
on the results of the High-z MareNostrum Simulation. From the numerical results 
we are able to estimate the observed restframe UV and Lymana emission in these 
galaxies. 

We find that our results are broadly consistent with the current observational 
constraints. The largest discrepancy is an overabundance in the faint end LF 
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MareNostrum High-z Simulation 
z=6, v=353 GHz 
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Figure 1. Flux in the ALMA band at 353 GHz as a function of the 
observed Lymana emission. The horizontal line represents the thresh- 
old for detection. In our model the brightest Lymana galaxies won't 
be detected by ALMA. 



both for LBGs and LAEs. The agreement between our model and observations 
for objects brighter than Mtjv < — 20 is remarkable. 

We also predict the expected flux at 353 GHz as a function of the observed 
Lyman-a emission. In our model, the brightest LAE galaxies are too faint FIR 
emitters. Correspondingly, bright IR sources will be detected by ALMA, but 
they will not be detected as bright LAEs. 
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